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Abstract 

We investigate the twist three distribution function e(x) in light-front Hamiltonian perturbation 
theory. In light-front gauge, by eliminating the constrained field, we find a mass term, an intrinsic 
transverse momentum dependent term and a 'genuine twist three' quark-gluon interaction term in 
the operator. The equation of motion relation, moment relation and the sum rules are satisfied for a 
quark at one loop. We compare the results with other model calculations. 
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Introduction 

Twist three parton distribution e(x) has not been explored broadly in the literature so far. e(x) 
is spin independent and chiral odd. Therefore it is difficult to measure it experimentally as it 



can combine only with another chiral odd object. It was shown in [l] that e(x) contributes 
to the unpolarized Drell-Yan process at the level of twist four. e(x) enters together with the 
chirally odd Collin's fragmentation function in the azimuthal asymmetry in semi-inclusive 
deep inelastic scattering (DIS) of longitudinally polarized electrons off unpolarized nucleons. 
This asymmetry was measured 2|-|4|, however, apart from e(x), several other distribution and 
fragmentation functions also appear in this asymmetry. In [5|, another possibility of measuring 
e(x) is proposed through two hadron production in polarized semi-inclusive DIS; where e(x) 
is coupled to a two hadron fragmentation function. The first calculation of e(x) was done in 

fl HQ 

[1], in MIT bag model. Further it was calculated in chiral quark soliton model p, uj, spectator 



model [8[ and perturbative one loop model (9]. The Q 2 evolution of e(x) has been calculated 
in jlO]. A nice review of the general properties of e(x) and the various model calculations can 



be found in [11]. In chiral quark soliton model, a Six) singularity was observed in e(x), which 
was later found also in a perturbative calculation [9j. However, no such 5(x) term was found 
in the spectator model as well as in MIT bag model. It is to be noted that as x = cannot be 
experimentally achieved, no direct experimental verification of the 5(x) contribution is possible. 
Another interesting aspect of e(x) is that the first moment of the flavor singlet combination 
of e(x) is related to the pion nucleon sigma term a^N The flavor non-singlet part of e(x) 
satisfies a sum rule connecting it to the hadronic mass difference between neutron and proton 
The first moment of e(x) obeys the sum rule 

J <b e(x) = ^(P,S\ mm I P,S). (1) 
The second moment of e(x) obeys another sum rule 

^ dx x e ( x ) = _ Nq . (2) 

where m is the mass of the quark and M is the mass of the proton. N q is the number of quarks 
of a given flavor. These sum rules are not satisfied in the bag model as the QCD equation of 
motion is modified in the bag. They are not satisfied in the spectator model either. In QCD 



equation of motion method (11], the first sum rule is saturated by the 5(x) contribution only, 
whereas in chiral quark soliton model, only a part of the contribution comes from the 5(x) 
term. In e(x) has been calculated for a quark dressed with a gluon at one loop in QCD in 
light front gauge. Starting from the Feynman diagram in A + = gauge, k~ was integrated out 
to get to perform the calculation in light-front time-ordered method. A 5(x) term was found 
which was shown to be related to the k + = modes. 

The disagreement between different model calculations for e(x) makes it an interesting object 
for further study. In this work, we calculate e(x) for a quark at one loop in QCD using light front 
Hamiltonian perturbation theory in A + = gauge. Instead of using the Feynman diagrams, we 
expand the state in Fock space in terms of multiparton light-front wavef unctions. The advantage 
is that these wave functions are Lorentz boost invariant, so we can truncate the Fock space 
expansion to a few particle sector in a boost invariant way. The two particle light-front wave 
functions (LFWFs) can be calculated analytically for a quark at one loop using the light-front 
Hamiltonian. The partons are on-mass shell objects having non- vanishing transverse momenta. 
They can be called field theoretic partons. The distribution functions can be calculated at 
the scale Q 2 using the multiparton LFWFs. To 0(a s ) this is an exact calculation. Another 
interesting aspect is the nature of ultraviolet (UV) divergence in e(x). Being a twist three 
object, e(x) contains a 'bad' component of quark field. As a result, it is expected to have a 
different UV divergence as compared to the logarithmic divergence in twist two distributions 
fi(x). In the following, we present our calculation. 

e(x) for a quark at one loop 

The twist three distribution e(x) is given by []| : 

P + f dy- ± xP+ 



e{x) = WJ ^ s ' WW*') 1 p > s) - (3) 

We introduce the projection operators A 1 * 1 = |7o7± and project out the light-front good and 
bad components of the quark field : 

^(±) = (4) 
In light-front gauge, A + = 0, the 'bad' component, ip^ is constrained, and the equation of 
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constraint is given by 12] 

,l>-(y-) = -^-(ia 1 ■ d 1 - + ga 1 - ■ A ± + /3m) V> (+) (jT); 
io + 



where the operator 4+ is defined as [12 ] 



1 _ 1 [°° 

d+f( x ) = i J d v e ( x -y )f(y ) 

The antisymmetric step function is given by 

e (x~) = V \ — e* ws 

vr y a; 

P denotes the principal value. For the dynamical field The operator in Eq. 

as 

a = m^y) = ^ Ht (o)7oV> (+ V) + ^ (+)t (o)7o^ H (^) 

Using the constraint equation for this can be written as: 

O e = O m + O g + o k - 



where 



O m = rra// +)t (0) 



1 - + a 



id + id + 



(5) 



(6) 



(7) 



can be written 



(8) 



(9) 



(10) 



o k = ^wt(o) 



d+ 



'ir, 



o 9 = # (+)t (o) [m-^) + Gr)- 



(12) 



O m is the mass term, O k is the transverse momentum dependent term and O g is the explicit 
quark-gluon interaction term, also called the 'genuine twist three' term. 



For we use two component formalism 
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(13) 



The two component field £(y) has the Fock space expansion 



dk+d 2 ^ 



b\(k)e iky + d„ x (k)e 



-iky 



(14) 



with 





(15) 



For a dressed quark state of momentum P and helicity a: 



P,a) =&&+(P j( r) |0) 



ci,A2 



dk + d 2 Jc- L 



2(2tt) 3 A;+ 7 v/2(2tt) 3 ^ 



2(27r) 3 P+<J d (P- fci - fc 2 ) 



i 2 (P, cr | fci, oi; fc 2 , A 2 )6 t (/ci, 0-1)^(^2, A 2 ) I 0). 



(16) 



Here a) and are bare gluon and quark creation operators respectively and <p\ and 2 are 
the multiparton wave functions. They are the probability amplitudes to find one bare quark 
and one quark plus gluon inside the dressed quark state respectively. We introduce Jacobi 
momenta x iy qi 1 - such that £^ %i = 1 and J2i Qi 1 ' = 0- They are defined as 

Xi - p+ ' 

Also, we introduce the wave functions, 



Qi 



hi a Xj P 



(17) 



i>i = <f>i, 
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which are independent of the total transverse momentum P ± of the state and are boost invari- 
ant. The state is normalized as, 



(P, A' I P, A) = 2(2TrfP + 5 x ,y5(P + - P ,+ )(5 2 (P 1 - P 



(19) 



The two particle wave function depends on the helicities of the electron and photon. Using the 



eigenvalue equation for the light-cone Hamiltonian, this can be written as 



i>2* 1 ,\( X >V L ) 



x(l — x) 



1 
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«»,,t 



p a x; 



- 2- 



cr 

x 



+ ima ± - 



x 



(20) 



m is t 
state 



le bare mass of the quark, a\ = o~2, 02 = — 01 ■ i^x actually gives the normalization of the 
13]: 



01 



1 a 'r 



l + x 2 , Q 2 

ax- log—. 

1 — X \x l 



(21) 



to order a s . Here e is a small cutoff on x. We have taken the cutoff on the transverse momenta 
to be Q. This gives the large scale of the process. The above expression is derived using Eqs 
f fT9|) . ffTB"]) and (120]) . In the above expression, we have neglected subleading finite pieces, /i is a 
small scale separating hard and soft dynamics such that (q^) 2 > y 2 - For a dressed quark, we 
get, 



x e k (x,Q 2 ) = 0; 



(22) 



xe 9 (x,Q 2 ) = ^ C f log-[-6(l 



x) 



l + x 



(23) 



x 



6(1 -x) 



l-e 



dy 



1 + y 2 



(24) 



where x e^x, Q 2 ), x e g (x, Q 2 and x e m (x, Q 2 ) are contributions from Ok,O g and O m respectively. 
In the above, we have used the normalization condition Eq. (f2T|) . One can write 



l + x 2 

1 — X 



-<5(l-x) 



dy- = + -5(1 - x 

1 — y 1 — x+ 2 



(25) 



The plus prescription is defined in the usual way, that is dxj 



fix) 



Qdx f{x) - f{l) . For a 

JO 1— x 

dressed quark state, M is the renormalized mass of the quark. The bare mass m of the quark 
is given in terms of the renormalized mass 14j : 



i rl , 3a s Q 2 
m — M [ 1 Cf log- 
Air /i 2 



(26) 



Using this, we can write 



a. 



xe(x, Q ) = 5(1 — x) + —Cflog 



Q 



2 r 



2x 



1 



1 - x + 2 XS{1 ~ X \ 



(27) 



27T " fl 

In the above result, the divergence at x — > 1 gets canceled by the contribution from the 



normalization of the state and we get the plus prescription. The above result agrees with 



9]. 



Note that there is no 5(x) contribution in xe(x). However, as shown in jj)], S(x) contribution 
in e(x) is related to the k + = modes (zero modes) in light-front gauge. By choosing the 
prescription given by Eq. OH]) we are avoiding the zero modes. xe(x) is zero if we take the quark 
mass to be zero. 

Sum rules and equation of motion relation 

xe(x) can be related to the twist two unpolarized quark distribution through the equation 



of motion relation 



15|: 



m 



xe(x) = xe(x) + j^fi( x ) 
Where fi(x) is the twist two unpolarized distribution function : 



fi(x) 



dy~ 



e-^ p+ y-(p,s\mi + Hy')\P,s). 



in 



(28) 



(29) 



Note that Eq. (|28[) is unaffected by the presence of the gauge link in the definition of the 
parton distributions. In the above relation, we have suppressed the scale dependence. e(x) is 
the genuine twist three quark-gluon interaction part of e(x). For a dressed quark, 



a* 



f 1 (x,Q*) = 6(l-x) + -!-C f log 

Z7T 



l + x 



(30) 



So the equation of motion relation is satisfied with x e(x, Q 2 ) = ^ Cf logQ? ^x5(l — x) — 1 + 
The first moment of e m (x) and e g (x) are given by 



x 



[ e m (x, Q 2 ) dx = 1 + ^ C f log%r (-loge - 1); 



(31) 



e g (x, Q 2 ) dx 



^C f log^ (loge + I 



(32) 



Each part has divergence as x — > 0. However, their total contribution is free of divergence : 



e(x, Q 2 ) dx 



f 1 a s Q 2 1 

/ (e m (x, Q 2 ) + e g (x, Q 2 ) dx = 1 + — ^ C f log— - 

Je 27T ll A 1 



2M 



[P,S\${0)r/i{0)\P,S). 



(33) 



So the sum rule in Eq. (TjQ) is satisfied. However, our x region is limited by and 1 as this 
is physically allowed. Note that this implies that the sum rule is not saturated by a 5(x) 



7 



contribution, unlike what was concluded in 11] . The possibility of a divergence in this relation 
was observed in jlj. The second moment of e(x) becomes 



f 1 Us 

l x e(x, Q 2 )dx = 1 Ct log 

Jo 27T 



Q 2 3 _ m 



(34) 



with 



1 ,1 

dx x e(x, Q 2 ) = dx 
o Jo 



X 



5(1 - x) - 1 + x 



0; 



(35) 



The rhs of Eq. (134j) vanishes in the chiral limit m — 0. Note that this result agrees with {9], 
the S(x) present there in e(x) does not contribute to this sum rule. The relation for the n-th 
moment of e(x) defined by [e] n = dx x n_1 e(x) can be written as jit] ]: 



m 



Nn = [e]n + J^lflL-V 

The n-th moment is calculated using the expression above : 



p Jo 



2tt 



n— 1 



2X 1 r/i \ 



2yr / %2L f-^j 2 



J'=l 



where we have used J Q dx 



l — X - 



x n " 1 t~z. . = — Eii i - On the r/is, 



12 /•! 



— C f log-^r / dxx n ~ 2 -5(1 -x)-l + x 
2tt J L2 



— C f log— 



2tt 



1 



1 



ll 



71—1 71 



(36) 



(37) 



(38) 



12 /■! 



n-2 



1 + x 2 3 _ 3i 
+ + 2 ~ 2 



n-l 



27r /r L j n Ti—1 



(39) 



So the moment relation Eq. fl36l) is satisfied. 

To summarize, in this Letter, we investigate the twist three distribution function e(x, Q 2 ) 
for a massive quark at one loop using light-front Hamiltonian perturbation theory in light-front 



S 



gauge. By expressing the operator in terms of dynamical fields, we find three terms in the 
operator ; quark mass term , intrinsic transverse momentum dependent term and an explicit 
quark-gluon interaction term. The intrinsic transverse momentum dependent part does not 
give contribution. The equation of motion relation directly relates xe(x, Q 2 ) to the twist two 
unpolarized distribution function fi(x, Q 2 ). The mass of the quark plays a vital role here. The 
first moment relation for e(x, Q 2 ) is satisfied without a 5(x) term in the distribution function. 
Contribution from the mass term as well as the quark-gluon interaction term diverge as x — > 0, 
however their total contribution does not diverge in this limit. The second moment vanishes in 
the chiral limit. The logarithmic divergence of e(x, Q 2 ) gives the scale dependence. 
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